General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



NASA Technical MemcNrarKium 73193 


(HHSA-TH*73193) SYSTEM OPTIHIZJITIOB OF M78-29036 

GASDfMAHIC LASERS, COHPOTER PROGRAM OSER'S 
HAHOAL (RASA) 49 p HC A03/MP A01 CSCL 20E 

Unclas 

G3/36 27236 


System Optimization of 
Gasdynamic Lasers 
Computer Program User's Manual 

L J. Often III, R. C. Saunders III, and S. J. Morris 


April 1978 


IVIASA 

National Aeror, tics and 
Space Admini, jfion 



NASA T^hnical Memorandum 73193 


System Optimization of 
Gasdynamic Lasers 
Computer Program User's Manuai 

L J. Otten III 

R. C. Saunders III, Ames Research Center, Moffett Field, California 

S. J. Morris, Langley Research Center, Hampton. Virginia 


NASA 

N Hional Aeronavilics and 
Sptice Administration 

Amm Research Center 

Moffut! F mid G<»iiforni(j n-tOdii 




SYSTEM OPTIMIZATION OF GASDYNAMIC LASERS 


COMPUTER PROGRAM USER' S MANUAL 

L. J. Otten III * and R. C. Saunders III** 
Ai^s Research Center 

and 

S. J. Morris 
Langley Research Center 

SUMMARY 


This report contains the user's manual for the computer program described 
in reference 1. Gasdynamic laser systems or subsystems can be optimized by 
using this manual in conjunction with reference 1. This manual assumes 
familiarity with gasdynamic laser systems and with large-scale computer models. 


PROGRAM DESCRIPTION 


Necessarily, the computer code that performs this analysis is large and 
complex. In its present form, it has 32 subroutines and in excess of 
6600 FORTRAN statements, of which perhaps 10-20% are comments. The code is 
written in a dialect of FORTRAN and was originally intended for the Control 
Data Corporation 6600/7600 computers. Hence, it uses a few language features 
and programming techniques that are unique to that machine. It has recently 
been converted to other machines, however, including the International Business 
Machines 360/370 computers, so its machine-dependence is minimal. A single 
analysis pass requires from 3 to 5 sec on the CDC 7600, and a typical optimiza- 
tion may require 300 sec '^r more. The program requires 106,000g words of 
memory on the CDC 7600 and operates without overlay or segments on this 
system. No external storage such as magnetic tape is required. 

The program consists of five modules, as shown in figure 1. The control 
module performs certain input and output functions, maintains the flow of 
information between modules, and calls the other modules at the appropriate 
time. The optimization module (CONMIN) is used to minimize or maximize some 
features of the system. The remaining three modules (GAIN, FLAME, and WTVOL) 
are the physical models and calculate the laser performance, combustion infor- 
mation, and weight/volume results, respectively. The modules and the FORTRAN 
subroutines associated with each of them are summarized in table 1. 


*Present address; 
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Air Force Weapons Laboratory, Kirtland Air Force Base, 
Wright Patterson Air Force Base, Dayton, Ohio. 
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TABLE 1.- 

PR(X!RAM MODULES AND 

SUBROUTINES 

Ifodule 

Function 

Associated 

subroutines 

CONTROL 

I/O control; 

C0PES2* 


conmunlcatlon 

COPE01 

COPE02 

COPE03 

COPE04 

COPE05 

COPE06 

COPE07 

CONMIN 

Optimization 

CONMIN* 

CNMN01 

CNMN02 

0^03 

CNMN04 

CNMN05 

CNMN06 

CNMN07 

CNMN08 

GAIN 

Laser gasdynamics 



and kinetics; gain, 

CYCLE* 


available energy 

CHAR 

DERFC 

FLUX 

TAUSUB 

RUNGEX 

EV 

NOZZLE 

CONTOUR 

MACH 

GAMMA 

TABNOZ 

AVG 

FLAME 

Combustion products 
flame temperature 

; FLAME* 

WTVOL 

Weight; volume 

WTVOL* 


^Denotes entry point. 
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Because this progrsn wss developed over an extended period of tlee, each 
of the physical aodel iMJdules (GAIN, FLAME, and WT\niL) was initially written 
to run on its own. Therefore, each of these modules was required to read the 
information It needed, execute, and write the output on Its own. For sim- 
plicity, this structure was retained in the integrated program. In addition 
to making module develofanent easier, the use of this structure permits easy 
substitution of any of the modules. Execution sequence is controlled by the 
CONTROL module. 

CoMaunicatlon between the modules and especially with the control module 
is handled by a single labeled COMMON block, GLOBCM. A parameter’s position 
in this block is used to identify its role in the sensitivity analysis or 
optimisation process. The contents of this block, called the global COhBION 
block, are shown In table 2. In the version described here, only the first 
3S elements of this block are used; the remaining 1465 are available for 
future use or can be deleted if memory space Is at a premium. 


TABLE 2.- GLOBAL CATALOG 


Location 

FORTRAN 

name 

Definition 

Units 

Defined by^ 

Used 

by 

1 

TSTAG 

Condiustor stagnation tempera- 

K 

1 . 

2 

1 . 

3 



ture 






2 

PSTAC 

Combustor stagnation pressure 

atm 


1 

1 . 

3 

3 

THROAT 

Cavity nozEle throat height 

mm 


1 

1 , 

3 

4 

AEXIT 

Cavity nozzle area ratio 

— 




1 

5 

XC 

Mole fraction of GO^ 

— 

1 , 



1 

6 

m2 

Mole fraction of H 26 

— 

1 . 

2 


1 

7 

XCO 

Mole fraction of CO 



1 


1 

8 

XHE 

Mole fraction of Hg 

— 


1 


1 

9 

ELSUB 

Length of subsonic portion 

cm 


1 


1 



of nozzle 






10 

ELNOZ 

Length of supersonic portion 

cm 


1 


1 



of nozzle 

cm 


1 



11 

C 

Length of 1st section of 

cm 


1 


1 



nozzle 






12 

RAO 

Length of 2nd section of 

cm 


1 


1 



nozzle 






13 

THETAl 

1 st nozzle angle 

deg 


1 


1 

14 

THETA2 

2 nd nozzle angle 

deg 


1 


1 

15 

BIG 

Maximum peak gain 

%/cra 


1 


1 

16 

PAVG 

Average available specific 

kJ/lb 


1 

1 , 

3 



energy 






17 

CHI 

Ratio of moles of carbon to 

— 


2 


■1 


moles of C-N-O-H fuel 


■^l - GAIN 

2 - FLAME 

3 - WTVOL 
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TABLE 2.- Concluded 


Location 

FORTRAH 

nwna 

Dafinltion 

Units 

Defined by^ 

Used by^ 

18 

Q 

Moles of excess oxidlser 

WWW 

2 


19 

R 

Ratio of moles of to moles 

of C-H-O-H fuel (Ol^AIR) 

WWW 

2 

2 

20 

S 

Ratio of moles of excess CO^ 
to moles of C-N-Q-H fuel 

WWW 

2 

2 

21 

DELKX 

Additional oxidizer temp, 
due to preheat 

K 

2 


22 

DELN2 

Additional makeup H., t^p. 
due to preheat 

K 

2 

2 

23 

QFRAC 

Mole fraction of Ng that Is 
preheated 

WWW 

2 

«> 

24 

WTO 

System weight (old bench) 

lb 

3 

3 

25 

VT 

System volume (old bench) 

ft^ 

3 

3 

26 

WT02 

System weight (new bench) 

lb 

3 

3 

27 

VT2 

System volume (new bench) 

ft* 

3 

3 

28 

B 

Cavlty/mlrror height 

in. 

3 

1. 3 

29 

ARCAV 

Cavity aspect ratio 

— 

3 

3 

30 

P 

Static pressure at end of 
cavity 

atm 

1 

3 

31 

EM 

Mach number at end of 
cavity 


1 

3 

32 

PC.AM 

Cavity process > 

— 

1 

3 

33 

ETA 

Diffuser efficiency (times 
normal shock) 

«ww 

3 

3 

34 

P0-P0WN 

Actual minus required pres- 
sure at end of dtffu.ser 
(>0 for recovery) 

psi 

3 

3 

35 

ETAN2 

Ratio of moles of Ng to 
moles of N, In air 
(OXMIR) 



** 


PROGRAM USAGE 


The program input Is divided Into five sections. The first section 
Is CONTROL and CONMIN input. The second section selects which options among 
the physical models will be used. The remaining sections .ire the inputs for 
the physical model raodtiles. The first two sections of Input are mandatory 
for any problem. The remaining sections of input depend on the p,irtlcular 
problem under investigation. 
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This section provides a card-by-card description of the Input and 
presents a sample. The format for each card is given In parenthesis after 
the contents* and a description of the contents Is given below it. 


A. CONTROL and CONMIN Input (Mandatory) 


Card 1: TITLE (20 a4) 

TITLE - any alphanumeric title* up to 80 characters. 

Card 2: NCALC, NDV, NSV, N2VAR, IPNPUT, IPSENS, IP2VAR (715) 

NCALC - calculation control. 

0: read input and stop. All other control input is optional. 

1: cycle through program once. All other control input is 

optional. 

2 1 do optimization. Additional control input of Cards 3-10 
is required. All other control input is optional. 

3: do sensitivity analysis. Additicnal control input of 

Cards 11-14 is required. All other control input is 
optional. 

4: do optimum sensitivity. Additional control input of 

Cards 3-14 is required. All other control input is 
optional. 

5: do two-variable analysis (all combinations of two variables). 

Additional control input of Cards 15-18 is required. All 
other control input is optional. 

6; do optimum two-variable analysis. Additional control input 
of Cards 3-10 and 15-18 is required. All other control 
input is optional. 

NDV - number of design variables in optimization or optimum 
sensitivity. 

NSV - number of sensitivity variables. 

N2VAR - number of objective functions in a two-variable analysis. 

IPNPUT - input print control. 

0: print card Images plus formatted print. 

1: formatted print only. 

2: no print of input. 



Card 3 
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IFSl^ - s«n«ltlvtt^ print control. 

Oi no detnilod print ««ch sensitivity anslysis. 
It detailed print each sensitivity analysis. 
(DEFAULT; 0) 

1P2VAE - two-variable analysis print control. 

0; no detailed print each two-variable analysis. 
1; detailed print each two-variable analysis, 

(DEFAULT; 0) 


orriMimio.s 

IPRIMT. ITMAX. tCNDIR, NSCAL, ITEM, LINOBJ, NACMX1» NFDG (815) 
IPRINT - optinleatlon print control, 

0; no print, 

1; print Initial and final opt Imitation Information, 

2; I plus function v.ilue and design variable values at e.ach 
Iteration. 

3; 2 plus constraint values, direction vector, and move 

parameter at each iteration. 

4; 3 plus gradient information, 

ITMAK - maximum numher of optimisation iterations (DEFAULT • 10) 

ICRDIR - conjugate direction restart parameter (DEFAIW - 3) 

NSCAL - scallttg frequency. Design v,ari,^bles are scaled every 
NSCAL iteration. (DEFAUIT - l», no scaling) 

ITRM - number of successive iterations which must satisfy conver- 
gence criteria (see Card 5) before optimisation is 
termiiuated. 

(DEFAULT - 3) 

LINOBJ - linear objective function Identifier. 

Oi objective function is nonlinear function of design variables, 
i: objective function Is linear function of design variables. 

(DEFAULT - 0) 


NACHXl - nufl^er of anticipated active constraints plus 1. 


(DEFAULT: NDV + 2) 

NFDG - finite difference gradient identifier. 

0: gradients computed by finite difference. 

1: gradient of objective computed analytically, gradients of 

constraints computed by finite difference. 

2: all gradients computed analytically. 

(DEFAULT - 0) 

Card 4: FDCH, FDCHM, CT, CTMIN, CTL, CTLMIN, THETA, PHI (8F10.2) 

FDCH - relative finite difference size for gradient calculation. 
(DEFAULT - 0.01) 

FDCHM - absolute finite difference size for gradient calculation. 
(DEFAULT - 0.001) 

CT - constraint thickness parameter. 

(DEFAULT - 0.05) 

CTMIN - minimum absolute value of CT 
(DEFAULT - 0.004) 

CTL - constraint thickness parameter for linear and side constraints. 
(DEFAULT “ 0.01) 

CTLMIN - minimum absolute value of CTL. 

THETA - mean value of push-off factor. 

(DEFAULT - 1.0) 

PHI - participation coefficient 
(DEFAULT - 5.0) 

Card 5; DELFUN, DABFUN (2F10.2) 

DELFUN - minimum relative change in objective function. Used to 
terminate optimization. 

(DEFAULT * 0.001) 
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DABFUN - nlnlnum absolute change In objective function. Used tc 
teralnate optlalaatlon. 

(DEFAULT - 0.001) 

Card 6; NDVTOT, lOBJ, SGNOPT (215.F10.2) 

NDVTOT - total number of variables linked to the design variables. 

NDVTOT must be greater than or equal to NDV (see Card 2). 
Tl.ls option allows two or more parameters to be assigned 
to a single design variable. The value of each parameter 
is the value of the design variable times a multiplier, 
which may be different for each parameter. 

(DEFAULT; NDV) 

lOBJ “ global catalog location of objective function for 
optimization. 

SCWPT - optimization control, 

-1: minimize objective function. 

+1; maximize objective function 

Card 7: VLB(I), VUB(I), X(I) (3P10.2) 

VLB(I) “ lower bound on 1th design variable. 

VUB(X) - upper bound on 1th oeslgn variable. 

X(I) - initial value of Ith design variable (If 4 0, It overrides 
the value read later) 

NOTE ; NDV (see Card 2) of these cards are required. 

Card 8; NDSGN(X), IDSGN(I). AMULT(I) (2I5.F10.2) 

NDSC^(I) - number associated » ith the Ith design variable. 

XDSC9j(l) - global catalog location of 1th design variable, 

AMULT(l) - constant multiplier on the Ith design variable. 

NOTE ; IIDV (see Card 2) of these cards are required. 

Card 9; NCONS ^^5) 

NCONS - number of constraint sets In an optimization. 
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’ J 

Card 10: IC0N(I)» JCON(I), LCON(I), BLU(I,1), BLU(I,2), BLU(I,3) 

(3I5.3F10.2) 

ICON(I) - first global catalog location of the Ith constraint 
variable. 

JCON(I) - last global catalog location of the 1th constraint 
variable. 

(DEFAULT - ICON(D) 

i.C0N(I) - constraint identifier for the Ith constraint variable. 

1: linear constraint 

0: nonlinear constraint 

BLU(I.l) - lower bound on Ith constraint variable. 

BLU(I,2) - upper bound on Ith constraint variable. 

BLU(I,3) - scaling factor on Ith constraint variable. 

(DEFAULT - (BLU(I.l) + BLU(I ,2))/2.) 

NOTE ; NCONS (see Card 9) of these cards are required. 


SENSITIVITY ANALYSIS 


Card 11: NSOBJ (15) 

NSOBJ - number of objective functions for sensitivity ar ilysls 

Card 12: NSENSZ(I). 1=1. NSOBJ (1615) 

NSENSZ(I) - global catalog location of ith sensitivity objective 
function. 


Card 13; ISENS(I), NSENS(l) (215) 

ISENS(I) - global catalog location of 1th sensitivity variable. 
NSENS(I) - nuari>er of values of 1th sensitivity variable. 

Card 14: SENS (I, J) 

SENS (I, J) - values of Ith sensitivity variable (J ■ 1 is 
nominal value) 


NOTE; NSV (see Card 2) of Cards 13-14 are required. 
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TWO-VARIABLE FUNCTION SPACE 

Card 15: N2VX. M2VX, N2VY, M2VY (415) 

N2)HC - global catalog 1. :atlcR of X-varlable. 

N2VX - number of valuea of X-varlable. 

N2VY - global catalog location of Y-varlable. 

K2VY - nm^er of values of Y-varlable. 

Card 16: N2VZ(I), I - 1, N2VAR (515) 

N2VZ - global catalog location of Ith objective. N2VAR is 
defined by Card 2. 

Card 17; XM2V(I), I - 1, M2VX (8F10.2) 

XM2V(I) - Ith value of X-varlable. M2VX Is defined by Card 15. 
viard 18: YM2V(I), I - 1, M2VY (8F10.2) 

YM?V(I) - Ith value of Y-varlable. M.’VY Is defined by Card 15. 
Card 19: EmJ 

The wrd ’END’ must appear here at the end of the control input. 

B. OPTION Input (Mandatory) 

A single card in NAMELIST format is used for inputting to this section. 
Card 20: WEIGHT, BURNER, LASER, VLCAV, DEBUG (5L, NAMELIST "OPTION") 

WEIGHT - Is the WTVOL module to be used? 

^DEFAULT - F) 

BURNER - is the FLAME module to be used? 

(DEFAULT - F) 

LASER - is the GAIN module to be used? 

(DEFAULT - T) 
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VTiCAV - la the cavity calculatloa to be carried 2 alrror dlamters 
downatream? If not, the Input value on Card 29 irlll be 
used. It la valid only If WEIGHT - T. 

(DEFAULT - F) 

DEBUG - write debug output? 

(DEFAULT - F) 

NOTE : The WTVOL module cannot be run alone (l.e., BUHNER =• F, 

LASER = F, WEIQIT - T), since certain parameters normally 
supplied by GAIN would not be available. FLAME and GAIN 
may be run alone. 


C. WTVOL Input 

A single card In NAMELIST format is used for inputting this data. 
This card should be present only if WEIGHT ■ T (see Card 20). 

Card 21: TP, TF, B, EMEXIT, P0EXIT, ETA, EFF, N, NOMISC, NOVALV (7F,1I,2L, 

NAMELIST "SIZE") 

TP - power out of laser cavity, In watts. 

(DEFAULT = 5 X 10 5) 

TF - shot length, in seconds. 

(DEFAULT * 2.0) 

B - nozzle/mirror height, in Inches. 

(DEFAULT -5.0) 

EMEXIT - Mach number at the diffuser exit. 

(DEFAULT - 0.9) 

P0EXIT - total pressure at the dlffuner exit, in PSI. 

(DEFAULT - 14.7) 

ETA - diffuser efficiency, in multiples of normal shock efficiency. 
(DEFAULT - 0.5) 

EFF - cavity extraction efficiency. 

(DEFAULT ’■1,0) 




(DKKAULT - 15, 



NOMISC “ In "roitiotfllani^om" wol|;ht Co b« omictod from Chi' total 
^MyMlom woitiht? 


UIKFAUI.T • V) 




NOVAIA’ - la "valve" weight to be oroiCft^d from the total ayatem ^ 

■O ' \ 

(OKKAltl.T - F) 


n. KIAMK Input 


A alngU card In NAMFI.TST format (a uaod lor Inputtiixg thin data. 
This card ahould be picaent only If lUiRNF.R ■ T (hoc Card 201, 


Card 22 ( 


CHI. Q. R, S, FTAN;, tJFRAC, KKUKL, KOXID 



(8F.21, NAMFI.IST "I'ROrKR") 


CHI - mo lea of C to mo Ion ui C-N-O-H. 
^ (PEFAia.T • 0.0) S_ 


Q - nh>le« of excraa air (11 all ia oxidlior). 

(OKFAUl.T - 0.0) ^ 


R - moles of N. to moles ol C-N-t>-H (It oxidtrer Is not air). .- 

D * 

S - mo lea of exoeua CO, to moles of C-N-O-H. 


(DKFAHl.T - 0.0) 


ETAN2 - moles ol to im'los of N, In air (If oxldlaer Is air). 



(OK FAULT - l.O) 

OFINZ - additional oxtdlier te*mperature due to preheat ing. In 
dogiees Kelvin. 

(OKFAlil.T - 0.0) 

0FI.N2 - additional nuikmip nttrogen temperature due to preheating. 
In degrees Kelvin. 


(hFFAm.T 


14 


QFRAC - fraction of makeup nitrogen at temperature DELN2. 

(DEFAULT • 1.0) 

KFUEL - fuel choice. 

1: carbon monoxide (CO) 

2; cyanogen (C .Nj) 

3: carbon and benzine (C + CgHg) 

4; carbon (C) 

5: benzonitrlle (C^H,^N) 

KOXID - oxidizer choice. 

1; oxygen (0,) 

2: nitrous oxide (N,0) 

3: nitrogen tetraoxide (N,0^) 

4; air 

E. GAIN Input 

This input section retains the format from the original code and, 
hence, a number of superfluous or redundant cards is still necessary, 
chiefly for (now nonexistent) plotting. Likewise, some of the input 
described below is for features of the code that have little utility in 
many problems, notably the restart capability. This entire section should be 
omitted if LASER * F (see Card 20). 

Card 23; LASTPLT ( 11 ) 

LASTPLT - an inactive plot parameter. It should be set to 1. 

Card 24; IDO, TITLE (II, 7A10) 

IDO - an inactive plot parameter. It should be set to 2. 

TITLE - any alphanumeric title, up to 70 characters. 

Card 25; TITLE2 

TITLE2 - any alphanumeric title, up to 80 characters. 

Card 26: TPLOT (TI) 

TPLOT - an Inactive plot parameter. It should be set to 3. 



IS 


III Vtmf (U) 

KfiST <* «n ln«cciv« flag. It ahould ba tat to 0. 

Card 28 1 COMITl (8A10) 

CMItl “ any alphanunaric title up to 80 charaetera. 

Card 29 i TSTAG» PSTAG, AEMT, THROAT, mX, AKCI^, FACTOR ( 7 F 10 . 0 ) 

TSTAG “ co«buator stagnation teii^eratur©, In dagroes Kalvin, 

If BlIRNKR •• T (sea Card 20 ) thia value la Ignored. 


PSTAG *• conbustor stagnation preaaure. In atmoapherea. 

AEXIT - nos*ie area ratio. 

THROAT - not, lie throat height. In wllllraeters. 

XrtAX « maxlmuw dowimtreaw distance over which cavity condlrjoos are 
computed, in centimetcra. If VLCAV " T and WEIGHT' • T 
(see Card 20), this value is ignored. 

ANGLE - cavity divergence angle, In degree*. 

FACTOR - nor.sl© friction factor (“1/2 pu^) 

Card 30l XC, J0I2, XCO, XHK, JJ (4F10.0,12) 

X€ mole fraction of CO,. If lUiRNER = T (sec Card 20), this 
value la Ignored. 

1012 — mole fraction of H .0, If lUIKNKR • T (see Card 20) and a 

watcr“^produC Ing fuel la iiBod , XH2 Is added to the reaulLlng 
water fraction. 


XCO - molt* fraction of He. 

dJ - rotational quantum number at which gain ts computed. 

NOTKi Gain Is also calculat ed at the retational quant unt minihet at 
which It is greatest. 


Card fl: 


ITAUSH, IKRTRf. iFl.liX, rSUg, ITAHU". IRATK, IH.’, Id. IPRINT 

(•>1 1 1 

JTAOSdt - water rate flag. 

0; uac nominal rate. 


1: use Stjarma late. 
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IRSTRT - restart flag. 

0; don't us® restart option. 

1: uae restart option. 

IFLUX - flux flag. 

0: cavity flux la not present. 

1: cavity flux Is present. 

ISUB - subsonic nozrle contour flag. 

0: use stored values. 

Ij use external contour (see Cards 40-41) 

2: us© CONTOUR subroutine (see Card 42) 

3: use previous values. 

4: special nosxle 1. 

5; special nozzle 2. 

6: special nozzle 3. 

7: use Mach number derivative subroutine, 

ITABLE - supersonic nozzle contour flag. 

0; use default contour (method-of-character 1st ics contour 
for a Laval nozzle) 

1: use table look-up. 

IRATE - vibrational probabilities flag. 

1: use external values (see Cards lA-lQl. 

0: use stored values. 

IH2 - hydrogen rates flag. 

1: don't use AFWL hydrogen rales. 

2; use AFWT. hydrogen rates. 

IJ - Inactive plot parameter. It should be set to 1). 

IPRINT - print (lag. 

0; normal print, 

1; additional print at every step I torn X = O.OX cm to 
X « 0.5 cm in the cavit' . 


■Sift 



o 

DELX - print Interval for cavity calculatlona, in cm. 

EI.SUB - length of converging '.subsonic) portion of the nozzle, 
in cm. A/A* - 10 at X - El.Sim. 

W ELNOZ - length of illverging (supersonic) portion of the nozzle . 

In cm. 

SLOPE - slope of the nozzle contour at X - 0. It is Ignored if 
ELNOZ - 0. 

O CURVE - radius of curvature of the nozzle contour at X • 0. 

Card 33: PHI. XBMON. XBMOFF (3F10.0) 

PHI - cavity optical flux. In watta/cm^. 

XBMi'N - downstream distance where flux starts. In cm. 


XBMOFF - downstream distance where flux stops. In cm. 

D 

NOTE : Include this card only If I FLUX ^ 0 (see Card 31). 

14: Bl. B2. B3. B4 



3St BS, B6. B7, B8 (4F18.0 

o 

17: m<) (K18.0) 

Bl-BH - coefficients of probability functions describing 

vibrational relaxation processes. They are specified 

B ln the form B(F| (T); where R( are constants and Fj 
are specific polynomials contained In subroutine TAUSl’H. 

NOTE: Include these cards only if IRATE - I (see Card 31). 

38: TAU31, SICC, SIGN, SIGN2 (IFIH.O 

39: SIGHK fF18.03 


TAU31 - transitional llfetiim' for the Cv'.iOOl) ► i;0,(1001 
tr.insltlon, in seconds. 


SICC - optical broadening collision cross section for collisions 
between CO,, in cm*. ** 


SIGN - optical broadening collision cross suction lor collisions 
between N,, in cm*. 


Card 

Card 


Card 


SsiGH2 - optical broadening collision cross section for collisions 
between HjO, .In cm^. 

SIGHE - optical broadening collision cross section for collisions 
between He, In cm^. 

NOTE ; Include these cards only if IRATE - 1 (see Card 31). 

40: CMl, CM2, CM3, CM4 (4F18.0) 

41; CMS (F18.0) 

CM1-CM5 - coefficients for polynomial fit to subsonic nozzle 
contour. 

NOTE ; Include these cards only If ISUB ■ 1 (see Card 31). 

42: C, R, THETAl, THETA2 (4F10.0) 

C - length, perpendicular to the throat plane, of the first 
section of the subsonic portion of the nozzle, in cm. 

R - radius of the circular arc defining the second section of 
the nozzle, in cm. 

THETAl - angle orienting the first nozzle section in degrees. 

The first section of the nozzle is at an angle of 
(90-THETAl) to the nozzle throat plane. 

THETA2 - angle orienting the second nozzle section in degrees. 

This section is a straight line, tangent to the circle 
at an angle of (90-THETA2) to the nozzle throat plane. 

NOTE ; Include these cards only if ISUB ■ 2 (see Card 31). 

43: N, X, EM, T (4F18.0) 

44: TVCl, TVC2, TVC3, TVN (4F18.0) 

45: TVH2, TVHH, P, A (4F18.0) 

N - mach step value at which calculation Is to be restarted. 

X - X coordinate at which calculation is to be restarted 
(corresponds to N) , in cm. 


EM - Mach number at X. 


T - gas temperature at X, in degrees Kelvin. 







TVCl - vlbratlonnl temperature at X of mode Vj of CO^, in degrees 
Kelvin. 


TVC2 - vibrational temperature at X of mode 
degrees Kelvin. 


TVC3 - vibrational temperature at X of mode 
degrees Kelvin. 


of CO^, in 


of COj, In 


TVN - vibrational temperature at X of N,, In degrees Kelvin. 

TVH2 - vibrational temperature at X of mode c*f H .0, In 

degrees Kelvin. 

TVHH - vibrational temperature at X of H 2 « in degrees Kelvin. 
P - static pressure at X, in atmospheres. 

A - nozzle to throat area ratio at X. 

NOTE! Include this card only if IRSTRT i 0 (see Card 31). 


This completes the input. Additional problems can be run by repeating 
Cards 1-41. To terminate, a card with 'STOP' is required. 



SAMPLE CALCinATION 


This section presents the results of a typical optlmi .-.ation problem. It 
utilizes virtually ail features of the model except that weight and volume 
results are not used. TItis calculation maximizes the available specific 
energy of a CO/NjO fuel /oxidizer device constrained to recover to atmospheric 
pressure. The design variables, constraints, and objective function are 
summarized in table 3. These cljoices are based on reference 2 and represent 
a reason^ibly advanced laser. 


TABLE 3.- OPTIMIZATION PARAMETERS 


Design variables 



Variable 


Initial Lower Upper 
value bound bound 


Stagnation pressure, atm 

60 

25 

115 

Nozzle throat height, mm 

.20 

0.15 

10 

Nozzle area ratio 

50 

10 

80 

Mole fraction of H,0 

.005 

.001 



Moles of Nj to moles of C-N-O-U 8 3.5 


15 












Constraints 


Constraint 


Stagnation temperature, K 

Peak gain, Xcm"^ 1*3 

Actual minus required pressure at end 314.5 
of diffuser (5^ Implies recovery) 


Ob.lectlve function 


Maximum available specific energy 



Shown below are the input data (fig. 2) and the resulting output (fig. 3) 
for this problem. 





GOL OPTIMIZATION SAMPLE, JULY 1 Q 77 
2 5 5 I 

5 05 3 T 3 0 

0 . 0 ? 0.01 

1 f> 1.0 

25.0 115.0 
0.15 «.^^ 10.0 

10.0 — 80.0 

0.001 0 . 0 *; 

3.5 15.0 

1 2 1.0 

2 3 1.0 ,, . 

3 ^ 1.0 

«♦ 6 1.0 

5 19 1.0 

3 

1 0 1000 . 2100.0 1000 . 

15 0 1.0 5.0 1 .0 

34 0 0 . 0 lOOO.C) lOCJl.O 

END 

SOPTIOK BURNFRsT, WETGHT = T, VLCAVsF, LASFt< = T, 

JSIZE ETA= 0 . 8 , EMEXIT= 0 . 2 , TEND 

fPPOP*‘R CHl= 0 ., 0 * 0 ., R= 8 ., S= 0 ., kFUEL* 


5 FNH 


1, KOXin=2, ^EN! 


20PTIHU7 GDL 
N2P OXICI/ER 
3 
0 

GOL OPTIMIZATION SAMPLE 
1 700.0 50.0 *;o.o 

0.13 0.005 0.0 

C005002C0 

2.0 0.1528 4.37 

STOP 


0.200 

0.0 


50.0 
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